Ultrafast light-induced magnetization dynamics in ferromagnetic semiconductors 
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We develop a theory of the magnetization dynamics triggered by ultrafast optical excitation of 
ferromagnetic semiconductors. We describe the effects of the strong carrier spin relaxation on the 
nonlinear optical response by using the Lindblad semigroup method. We demonstrate magnetization 
control during femtosecond timescales via the interplay between circularly polarized optical excita- 
tion, hole-spin damping, polarization dephasing, and the Mn-hole spin interactions. Our results 
show a light -induced magnetization precession and relaxation for the duration of the optical pulse. 
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With the discovery of III-Mn-V ferromagnetic semi- 
conductors 1], one can envision a new class of multi- 
functional devices that combine information processing 
and storage on a single chip with low power consump- 
tion. One of the challenges facing such future devices 
concerns the speed of the basic processing unit, governed 
by the dynamics of the collective spin and its ability to 
respond quickly to time-dependent external fields. Ultra- 
fast pump-probe magneto-optical spectroscopy has been 
used to study the magnetization dynamics triggered by 
femtosecond optical pulses in both metals 0,0 and ferro- 
magnetic semiconductors 0, 0, IE Q • Transient magnetic 
effects have been observed even in the initial nonthermal 
temporal regime, where the concept of carrier tempera- 
ture is not meaningful 0,0] • 

In III-Mn-V semiconductors, the ferromagnetic order 
is induced by the interaction between the itinerant va- 
lence band hole spins and the localized Mn spins 0, El • 
The magnetization is therefore sensitive to the itinerant 
carrier properties, which in III-V semiconductors can be 
well controlled with ultrashort optical pulses The 
combined magnetic and optical properties of this sys- 
tem open new possibilities for manipulating the collec- 
tive spin without using magnetic fields. This can be 
achieved, e.g., by exciting spin-polarized carriers with 
circularly-polarized femtosecond optical pulses. The ef- 
fects of such photoexcited spin on the ferromagnetic state 
can be studied by monitoring the change in the macro- 
scopic magnetization as function of time using ultrafast 
spectroscopy 0, 0, 0, 0, El • Light-induced magneti- 
zation and very fast Mn spin rotation towards the di- 
rection of propagation of the circularly polarized optical 
pulse, perpendicular to the ground state magnetization, 
were reported in Ga(Mn)As ferromagnetic epilayers Q. 
Kimel et.al. observed a photoinduced magnetization 
attributed to the photoexcited conduction electron spin. 
Kojima et.al. 0] interpeted their data in terms of a rapid 
(ps) rise in carrier temperature, followed by a slow rise in 
the spin temperature and a magnetization reduction over 
100s of ps. For very intense photoexcitation of In(Mn) As, 
Wang et.al. 0] observed a fast (100s of fs) and a slow 



(100s of ps) demagnetization regime. Other possibilities 
for light-induced magnetic effects in non-magnetic semi- 
conductors were recently explored theoretically 0, 0] . 
In doped semiconductors, a light-induced Kondo effect 
was predicted [l2T |. while a light-induced paramagnetic- 
to-ferromagnetic phase transition was predicted in un- 
doped semiconductors fl3| . 

In this letter, we develop a theory of the ultrafast 
nonlinear optical dynamics and light-induced hole and 
Mn spin relaxation in ferromagnetic semiconductors. 
We treat the magnetic exchange interactions within the 
time-dependent Hartree-Fock approximation, and treat 
the hole spin damping and resulting nonlinear polariza- 
tion dephasing within the Lindblad semigroup method 
fl~i|. We predict an ultrafast light-induced precession 
and relaxation of the collective magnetization towards a 
direction determined by the nonlinear optical polariza- 
tions. The above spin dynamics results from the inter- 
play between the nonlinear optical excitation, the inter- 
action between the photoexcited hole and Mn spins, and 
the strong hole-spin damping and optical polarization 
dephasing, The predicted dynamics is observable with 
ultrafast magneto-optical pump-probe spectroscopy. 

We adopt the simple 0, 0] two-band sp-d model 
Hamiltonian H(t) = K e + Kh + -ff e xch + Hr,{t), where Kh 
(K e ) is the kinetic energy of the valence band heavy holes 
(conduction band electrons) created by h\^ a (e^) with 
total angular momentum J z = ±3/2 (spin a = ±1/2) 
and dispersion e£ — k 2 /2rrih (e£ = k 2 /2m e ) (we set 
h=l). Unlike in II- VI semiconductors, the Mn impu- 
rities provide both a hole Fermi sea and randomly dis- 
tributed S — 5/2 localized spins, Sj, located at positions 
Rj . These local spins interact strongly with the itinerant 
hole spins via the Kondo-like exchange interaction 
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where a is the hole spin operator and V is the volume. 
The coupling of right-circularly polarized optical holds is 
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described in the rotating frame by 

H L (t) = -d(t)J2€iht^ + h.c, (2) 



where d(t) = fi cv £(t) is the Rabi energy, \x cv is the dipole 
transition matrix element, and £(t) oc exp [— i 2 / T p] is 
the optical pulse, with central frequency us p and duration 
T p . In the above Hamiltonian, the spin quantization z- 
axis was taken parallel to the direction of optical pulse 
propagation, which similar to Ref. |7j was chosen here 
perpendicular to the ground state spin polarization. 

The ground state of the above Hamiltonian is typically 
described by treating the interaction Eq.Q) within the 
mean field virtual crystal approximation, which neglects 
spatial correlations and assumes uniformly distributed 
classical Mn spins 0, 0|. Within this approximation, 
the holes experience an effective magnetic field propor- 
tional to the average Mn spin The valence states 
then split into two spin-polarized bands, separated by 
the magnetic exchange energy JcS, where c is the Mn 
concentration, and populated by the hole Fermi sea. The 
spins of the above bands point antiparallel and parallel 
to the Mn spin, which is taken here to define the x-axis. 

We calculated the ultrafast nonlinear optical response 
by solving the equations of motion for the density ma- 
trix. For right-circularly polarized light, spin-| electron 
states are not populated for the above Hamiltonian, and 
we only need to consider the interband optical polariza- 
tions P kcr = (/i-ko-eki), the carrier populations, the aver- 
age Mn spin S, and the hole spin s£, whose components 
are S L = E CT v(hi a h^) and s£ + = (hlyh kl ), where 
s k+ = s ka . + is^y. By factorizing all higher density ma- 
trices, we obtain a closed system of equations analogous 
to the semiconductor Bloch equations 01 generalized to 
include the spin degrees of freedom. Due to the selection 
rules, Eq. J2J), right-circularly polarized light excites di- 
rectly the spin-| hole-spin-! electron polarization Py. 
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where S±(t) = S x (t)±iS y (t) and ft k = el+e^-oj p -i/T 2l 
where Ti is the polarization dephasing time. The rhs 
of the above equation describes the Pauli-blocking non- 
linearity (Phase Space Filling) [ll|, determined by the 
equations of motion for the carrier populations, with ini- 
tial condition the Fermi-Dirac distribution of the hole 
Fermi sea. Additional nonlinearities come from the light- 
induced changes in the Mn spin S(t) from its ground state 
configuration, discussed below, which result in a time- 
dependent effective magnetic field JcS(i). The latter 
gives a time-dependent e-h pair energy oc S z (t), which 
renormalizes the spin-polarized hole bands, and a time - 



dependent coupling oc S-(t) between Pj and Pj 



t - O k + ^ W - ^P^ = d(i) S £ + . (4) 



Pi is also generated due to the hole spin coherence, 
photoexcited or ground state, described by si (rhs of 
Eq.Q). In turn, the hole and Mn spins depend on the 
optical polarizations and precess due to the exchange in- 
teraction as described by the equations of motion 

d t s£ = JcS x s£ + d(t) h k , d t S = JJ2 s k x s -(5) 



In the above equation, the photoexcited hole spin is de- 
termined by the nonlinear optical polarizations (and thus 
the interactions) as described by the vector 



h k = Im (P kl ,-iP ki ,P kT ). 



(6) 



During the femtosecond time scales of interest here, we 
neglected for simplicity the much slower precession and 
relaxation of the Mn spin due to the effective magnetic 
field caused by the magnetic anisotropy and the Gilbert 
damping [Tjilig. The nonlinear effects were treated non- 
perturbatively by solving numerically the above system 
of coupled equations. 

Relaxation in ferrmagnetic semiconductors mostly oc- 
curs during time scales longer than the pulse duration. 
A notable exception is the hole spin relaxation, which 
occurs within 10's of fs due to the strong spin-orbit 
coupling in the valence band and the disorder-induced 
scattering between the different momentum states [l7j . 
Carrier spin relaxation is typically described within the 
spin Bloch equations by introducing relaxation terms so 
that the spin relaxes to a direction opposite to the Mn- 
spin [id ll^ . However, this spin relaxation also leads 
to strong dephasing of the optical polarizations and cou- 
pling between Pp and Pj_. Thus, the dephasing of both 
the spin and the optical polarizations must be treated 
consistently, including the light-induced nonlinear con- 
tributions. We treated all the above effects by using the 
Lindblad semigroup description of dissipative quantum 
dynamics [l^j . Under the general assumptions of linear 
coupling between bath and system operators L k and posi- 
tivity and semigroup-type density matrix time evolution, 
the relaxation contribution is given by [l^j 
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The Lindblad operators L k = h^h^, where h J k ^ ( h^) 
creates a hole with spin parallel (anti-parallel) to the 
Mn spin, describe relaxation of the hole spin toward the 
direction opposite to S(t) with a rate r kk ' = S^Tg. 
The detailed equations, obtained from Eq.© by fac- 
torizing the higher density matrices, will be presented 



FIG. 1: Mn-spin trajectory for three different relaxation 
regimes: (a) T 2 =2Ti=330 fs, T s = 0, (b) T 2 =2Ti=330 fs, 
l/r s =10.5 fs, (c) 7W65 fs, T 2 =l/r s =10.5 fs. 



elsewhere. Here we note that the hole spin compo- 
nent perpendicular to S dephases with a rate T s , while 
the spin component parallel to S relaxes with a rate 
2r s towards -s max + s^ ax - (s£) 2 , where s max , with 
Smax — s max = n k/^ ~ ("-k/^) 2 , is the maximum hole 
spin for given number of holes n£ = J2a(^t.a^ a )- Eq-Q 
also gives optical polarization dephasing and coupling be- 
tween P-j- and P| that depend on S(t), s h (t), and the hole 
populations, whose expressions will be given elsewhere. 
We note that electron and Mn spin relaxation can also 
be included by using analogous Lindblad operators. 

Fig. 1 shows the light-induced time evolution of the 
Mn spin for pulse duration r p = 250fs, Rabi energy d 
= 60meV, and detuning 10 meV below the threshold of 
interband absorption. Other parameters used in our cal- 
culation are the exchange energy JcS= 125 meV, hole 
Fermi energy 100 meV, fraction of initial holes 0.33 of 
the Mn impurities, and m^/m e = 7.15. To elucidate the 
role of relaxation, we compare in Fig. 1 the Mn spin 
trajectories for different values of T s and T2. The Mn 
spin rotates in a clockwise direction within the x-y plane 
perpendicular to the direction of optical pulse propaga- 
tion. With increasing hole spin damping T s , it develops 
a large out of plane z-component that increases further 
with decreasing T 2 . 

The spin dynamics of Fig.l comes from the Mn spin 
precession around the mean hole spin. Unlike for a mag- 
netic field, the latter changes strongly with time, as de- 
termined by the nonlinear optical polarizations, the spin 
relaxation, and the precession around S (Eq.JSJ)). For 
the parameters in Fig.l, the hole spin follows overall the 
Mn spin trajectory, with the exception of a component 
perpendicular to S, determined by the nonlinear optical 
polarizations and the spin relaxation, which triggers the 
Mn spin dynamics. The light -induced deviations from 



FIG. 2: Trajectories of the light-induced field h (units of 
(1(0)/2Ef) for the three relaxation regimes of Fig. 1. 

the equilibrium spin configurations occur within the pulse 
duration and increase with its duration and intensity. 

At first glance, the clockwise Mn spin rotation is 
suprising given that, due to the selection rules Eq.(0), 
right-circularly polarized light creates spin-up holes that 
should lead to counter-clockwise rotation (see Eq.JSJ). In 
fact, our numerical results show that the z-component 
of the photoexcited hole spin is negative due to the ex- 
change interaction. This is a consequence of the fact that 
s h is photoexcited via the nonlinear optical polarizations 
and not by the optical field directly (see Eq.JSJ)). Also 
striking is the strong effect of the hole spin relaxation 
on the Mn spin trajectory: it leads to the development, 
within the pulse duration, of a large component S z in 
the direction of optical pulse propagation despite the ab- 
sence of any Mn spin damping in Eq. |JSJ . A very fast Mn 
spin tilt toward the z-direction was observed experimen- 
tally 0, however for weaker photoexcitation. The spin 
dynamics in Fig.l is significant for Rabi energies larger 
than a few tens of meVs, which corresponds to stronger 
photoexcitation than in Ref. [4 but not as strong as in 
Ref. 0|, where heating effects^} become pronounced. 

To interpret the above results, we consider the sta- 
tionary (adiabatic following) limit, where the hole spin 
precesses and/or relaxes much faster than the pulse du- 
ration and period of Mn spin precession, and can thus 
adjust to the instantaneous d(t) and S(i) values. Ne- 
glecting the component of d t s h perpendicular to S, we 
can then eliminate the hole spin from Eq. (J5J and obtain 
after some algebra and noting that S • dtS = 

dtS = fj 1 ® p2 [r s h x S - Jc S X (S X h)] , (8) 

where h(t) — ^ k is an effective magnetic field de- 
termined by the nonlinear optical polarizations (Eq.©). 
The above equation for the Mn spin demonstrates a light- 
induced precession around h(i) and Gilbert damping-like 
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relaxation of S toward h(t). Such ultrafast nonlinear 
optical effects should be contrasted to the static effects 
governed by the magnetic anisotropy (lol Ht| . 

The light-induced effective magnetic field h(i), Eq.©, 
that governs the Mn spin dynamics has z-component de- 
termined by the dissipative part of the nonlinear optical 
polarization, ImPj , and in-plane components determined 
by Pj_. Its time-dependence is shown in Fig. 2. The ex- 
change interaction plays an important role in determin- 
ing h(i). To see this, we note that, for J = 0, P; = 0. 
The corresponding h then points along the z-axis, with 
magnitude determined by the dephasing of Pf and, for 
T% < T p-, proportional to the optical pump rate. In this 
case, the Mn spin would relax to the z-axis. However, the 
exchange interaction and spin relaxation drastically af- 
fect the optical polarizations and thus h(i). For example, 
by coupling the hole spin-f and spin- J. states, they lead 
to a large P^ and thus a large x-y plane component of the 
effective magnetic field h. This is demonstrated in Fig. 2, 
which shows h(t) for different values of the relaxation 
rates. In fact, for sufficiently long T 2 and detunings close 
to or below the interband absorption threshold, ReP| 
exceeds ImPj , in which case the light-induced magnetic 
field h points mostly along the x-y plane, in a direction 
almost perpendicular to the ground state magnetization 
(x-axis) (see Fig. 2(a)). For T s = 0, the second term on 
the rhs of Eq. ||SJ then leads to Mn spin relaxation mostly 
within the x-y plane, while with increasing T s , the first 
term leads to precession out of the x-y plane, consistent 
with the numerical results of Fig.l. The polarization de- 
phasing, due to the hole spin damping as well as the 
disorder and other contributions, enhances the relative 
magnitude of ImPf and thus h z (see Figs. 2(b) and 2(c)). 
As a result of h z , S develops an additional z-component 
as it relaxes toward h. One should note that the rhs of 
Eq. JEJ vanishes after the optical pulse is gone and thus 
the light-induced Mn spin relaxation and precession only 
occurs during the pulse. Furthermore, the optical polar- 
izations, and therefore h(i), depend on the spins, and 
thus the Mn spin dynamics leads to an additional time- 
dependence of h(t) and to corresponding nonlincarities. 
This can be seen in Fig. 2, where the trajectories of h(t) 
have a loop-like shape. Finally, the Phase Space Filling 
optical nonlinearity also affects the time-dependence of 
h(t), Eq.©, and thus the Mn spin relaxation Eq.©. 

To gain insight into the effects of the optical excitation 
and exchange interaction on the hole spin dynamics, one 
can solve analytically the polarization equations of mo- 
tion, Eqs.Q and Q, for < t p in the stationary limit, 
discussed elsewhere. For example, we obtain this way for 
the x-y spin components (r s = 0) 
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where = ^ a $k a hv.o) ■ The first term on the rhs de- 
scribes hole spin precession around S(t) triggered by a 
time-dependent effective exchange interaction. The sec- 
ond term describes a light-induced in-plane rotation of 
the hole spin, even for S z — 0, while its dissipative part 
describes a light-induced hole spin relaxation. Finally, 
we note that, for very strong photoexcitation [j| and/or 
detunings well above the absorption edge, the heating of 
the hole Fermi sea can become significant. Similar to 
metals Q, such hot hole effects can be treated by intro- 
ducing a time-dependent effective temperature. 

In conclusion, we presented a theory of ultrafast non- 
linear optical and spin dynamics in III(Mn)V ferromag- 
netic semiconductors. We demonstrated a light-induced 
Mn spin precession and relaxation, which can lead to a 
magnetization tilt of several tens of degrees. We showed 
that the latter is determined by the nonlinear optical po- 
larizations and depends critically on the exchange inter- 
actions, dephasing, and hole spin relaxation. The non- 
linear light-induced magnetic effects discussed here occur 
during the optical pulse and can be observed with ultra- 
fast pump-probe magneto-optical spectroscopy. 

This work was supported by the EU Research Training 
Network HYTEC (HPRN-CT-2002-00315). 
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